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Abstract: Odonata (dragonflies and damselflies) are efficacious for management and conservation 
efforts in freshwaters. In recent times, increased effort has gone into enhancing awareness, data and 
information on dragonflies among scientists and policymakers. Here, we examined the Odonata 
community of dams within the Matobo National Park, a world heritage site in southern Zimbabwe. 
Specifically, we determined diversity of larval and adult odonate in five dams over one year. Several 
physico-chemical parameters were measured, including pH, total dissolved salts (TDS), electrical 
conductivity (EC). Further, we assessed how habitat types (i.e., detritus, sand, gravel, plant type) af- 
fect Odonata communities. In addition, surveys were conducted to identify critical threats to Odo- 
nata in the Matobo National Park and surrounding areas. Broadly, results showed that Odonata 
nymphs and volant adults did not change predictably across all sampling occasions. Further, pH was 
positively correlated with Libellulidae and Platycnemididae, whereas TDS and EC were positively 
correlated with Libellulidae, Platycnemididae and Macromiidae. Contrariwise, Gomphidae were 
found to be negatively correlated to TDS, EC and pH. The threat analyses identified natural system 
modifications, agricultural expansion and intensification, as well as human intrusions and distur- 
bance as the major threats to Odonata and freshwater resources in Matobo National Park. Taken to- 
gether, these data provide baseline data that will be useful for future monitoring of threats and sub- 
sequently conservation strategies in the Matobo National Park and other protected areas in Southern 
Africa. 
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Introduction 


Odonates (dragonflies and damselflies), as one of the earliest winged insects from the middle Perm- 
ian (Kalkman et al., 2008) are prominent in many types of freshwater habitats globally (Dijkstra et 
al., 2014). Owing to their complex life histories, odonates are good bio-indicators of the overall health 
of an ecosystem in both terrestrial and aquatic habitats. Since they inhabit both aquatic and terres- 
trial environments, they may better reflect environmental variations in different ecosystems. 


Dragonflies and damselflies have been successfully used as major components in assessing the eco- 
logical condition of aquatic biotopes in some African countries (e.g., South Africa and Kenya; Claus- 
nitzer, 2003; Suh & Samways, 2005). Dragonflies and damselflies have also been promoted as an 
easy-to-learn tool in environmental education programs. For instance, the Dragonfly Biotic Index 
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(DBI), a biomonitoring tool that uses the presence of dragonfly at a site to define the health of that 
freshwater ecosystem, was developed for South African freshwater systems (Simaika & Samways, 
2009; Samways & Simaika, 2016). The aforementioned works have paved the way for projects on 
water conservation using bioindicators, especially in developing countries, where pressures on fresh- 
water systems are prevalent (Clausnitzer et al., 2017). 


The recognition that Odonata species as useful indicators of ecosystem health in environmental im- 
pact assessments and monitoring programs, has seen an increased effort in making information on 
dragonflies available to scientists, conservation area managers and policymakers (Cardoso et al., 
2020; Cordoba-Aguilar, 2008). However, effective use of odonates in conservation efforts, requires 
baseline information on the diversity and abundance of these species that are currently understudied 
in most of the African continent (Rodriguez-Tapia et al., 2020; Samways & Steytler, 1996). 


The distribution of various species of Odonata is highly variable across spatio-temporal scales. In 
Africa, 15 families have been recorded (Dijkstra & Clausnitzer, 2014; Dijkstra et al., 2015) with evi- 
dence that some families are restricted to cool streams or rivers, while others are restricted to ponds 
or still clear waters, and some to marshy places. High diversity of Odonata species is often at sites 
that offer a wide variety of microhabitats. 


Many biotic and abiotic factors affect the distribution of aquatic insect larvae (sensu Schriever, 
Cadotte, & Williams, 2014). For instance, acidity of the water, the amount and type of aquatic vegeta- 
tion, the temperature, and water flow all affect the distribution of Odonata larvae and can all poten- 
tially affect the distribution of aquatic insect larvae (Schriever, Cadotte, & Williams, 2014). While 
some species are affected by a myriad of biotic and abiotic factors, some species can tolerate a broad 
range of conditions (Acquah-Lamptey, Kyerematen, & Owusu, 2013). Broadly, dragonflies tend to be 
much more sensitive to pollution than damselflies (Acquah-Lamptey, Kyerematen, & Owusu, 2013). 


Dragonflies occupy a wide range of habitats with many species and families having their specific 
environmental preferences. The habitat association of dragonflies and damselflies in Africa has been 
described by many researchers (e.g., Deacon, Samways, & Pryke, 2020; Vorster et al., 2020). More- 
over, several studies have shown that species assemblages of odonates are influenced mainly by 
aquatic and terrestrial vegetation (Paice, Chambers, & Robson, 2017). Studies have shown that adult 
diversity and abundance is often positively correlated with local abundance of vegetation (Remsburg 
& Turner, 2009). The adults of some species require emergent or waterside plants to use as perches 
(Pritchard & Kortello, 1997; Worthen & Morrow, 2016). For instance, McCrae & Corbet, (1982) 
documented that the adult males of Tetrathemis polleni (Selys) will compete for the tops of the sticks 
as perches will keep returning to the same perch even after a disturbance. Other species such as The 
Large Red Damsel (Pyrrhosoma nymphula) need specific submerged or floating macrophytes to lay 
their eggs on (Rehfeldt, 1990). Different odonate species can also occupy the same small area if each 
species can use a different substrate for perching and oviposition (Martens, 2002). 


To date, there are many taxonomic studies on odonate larvae and imagoes of Zimbabwe (Pinhey 
1961; Pinhey 1984; Tarboton & Tarboton, 2019) and an extensive collection of odonatan species col- 
lected around Zimbabwe by world renowned collectors like Dr. Elliot Pinhey, whose samples are all 
deposited in the Natural History Museum of Zimbabwe. Some of the aforementioned collections 
carry useful occurrence records that have been used in predictive modelling for the distribution of 
endemic species such as Chlorestes elegans (unpublished data, Mafuwe 2020). Nonetheless, most 
records date back at least 30* years (Pinhey 1961; Pinhey 1984; Kipping, Dijkstra, Clausnitzer, Suh- 
ling, & Schütte, 2009; Dijkstra & Clausnitzer, 2014). A few studies in Zimbabwe have also focussed 
on quantitative habitat correlates of diversity and habitat associations. For example, researchers have 
assessed the effect of industrial waste and temperature fluctuations on Odonata diversity in Zimba- 
bwe (Dube, Makaka, & Sibanda, 2010; Makaure et al. 2015). However, the ecological studies on 
richness and distribution at national levels remain few; these studies are also few in lacustrine sys- 
tems, especially those in protected areas. 


Strikingly, in Zimbabwe as in many other parts of sub-Saharan Africa, biodiversity conservation 
priority lists, omit Odonata (Burgess et al., 2005; Mutasa & Ndebele-Murisa, 2015), even though they 
are known to occur in these conservation areas (Bernard & Daraz, 2018; Dijkstra & Clausnitzer, 
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2014). Considering, a few species are listed as threatened or data deficient in the IUCN Red list 
(IUCN, 2017) cursory data on odonate diversity and distribution will be useful for the conservation 
of odonate species, along with their aquatic habitats. 


Determining the factors that may affect the presence of a species is crucial to determine the species 
that can be affected by environmental or habitat changes (Cardoso et al., 2020; Nagy et al., 2019). In 
this study, we describe shifts in adult and larval Odonata abundances in an understudied world heri- 
tage site in Zimbabwe over one year (across three seasons), as well as determine the relationships that 
exists between water quality, aquatic habitats and these odonate fauna. Specifically, we aimed to (i) 
describe the odonate fauna in terms of abundances and diversity over a one-year period, (ii) assess 
relationships between Odonata and key environmental and habitat variables and, (iii) assess current 
and future threats to Odonata diversity. Along with these objectives we predicted that diversity of 
odonate nymphs would be higher in the cold seasons where emergence of odonates is low. Con- 
versely, we predicted that volant dragonflies would have higher diversity in the warmer months when 
emergence of odonates is high. Outcomes from this research are expected to form the basis for Odo- 
nata research and subsequently contribute to conservation actions through long-term ecosystem 
monitoring of aquatic fauna. This study is also a first step towards identifying bioindicators urgently 
needed for wetland monitoring in Zimbabwe. 


Material and methods 


Study Area 


The Matobo Hills in Zimbabwe was declared a World Heritage Site in August 2003 and is character- 
ized by a rocky landscape with a unique natural and social heritage (Sagonda & Pegg, 2015). The 
Matobo Hills comprise of wooded valleys and boulders and are a place of very high biodiversity of 
flora and fauna (Matobo Hills Management Plan, 2015). The world heritage site covers an area of 
about 3100 km?, of which 424 km’ is a National Park (Figure 1), the remainder being mostly com- 
munal land and a small proportion of commercial farmland (Gambiza & Nyama, 2006). Part of the 
national park is set aside as 100 km? of game park. The park also extends along the Thuli, Mtshelele, 
Maleme and Mpopoma river valleys (Figure 1). 


Sample collections 


We collected samples for Odonata nymphs and adults at the start, mid and end of each season across 
five dams (Figure 1) between May 2016 and April 2017. Sample collection was designed to coincide 
with Zimbabwe’s three climatic seasons (hot wet season from mid- November to March, a cold dry 
season from April to July, and a hot dry season from August to mid-November). Sampling was con- 
ducted between 10:00 hours and 16:00 hours when adult Odonata are mostly active. 


At each of these dams, six sampling points were selected within 100 meters of the main access points 
and were at least 200 m apart. Using a South African Scoring System (SASS) net (35 by 35 cm and 
mesh size of 0.8 mm), odonates were disturbed and encouraged to drift by agitating the substrate 
using the traveling kick method for a 5-min period over a 5 m2 area of the pool. Contents from the 
aquatic nets were emptied into a pan to pick out the nymphs using forceps and a total of six scoops 
were done at each sampling point. Specimens were preserved in ethanol | and transported to the 
laboratory for counting and identification to families and morphospecies. At each point recordings 
were made on micro-habitats according to the presence or absence of aquatic plants as well the bot- 
tom substrates (i.e., detritus, sand, gravel) from which the nymphs were collected. The following 
water quality parameters were measured at each point using a multiparameter meter (YSI Inc, USA): 
pH, Total dissolved salts, electrical conductivity. 
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Figure 1. Study area (Matobo National Park) showing five dams sampled across three seasons (hot wet, 
cold dry, hot dry) between May 2016 and April 2017. 


For observations of adult specimens, sampling was done along two, 500 m transects on both sides 
of the dam, with a width of 5 m perpendicular from the water systems. Odonata species encountered 
were identified and counted. Identification of species was primarily done using close-focusing bin- 
oculars and field guides. Where necessary one individual of each species was caught and examined 
using a hand-lens to confirm identification and subsequently released. Species were photographed 
to provide a permanent record of identification. Where further identification was required, voucher 
specimens were captured and transported in labelled paper envelopes for further identification us- 
ing keys by Dijkstra & Clausnitzer (2014). The Dragonfly Biotic Index (DBI) for each species was 
also noted. The Dragonfly Biotic Index is a compound index based on three criteria: geographical 
distribution, conservation status and sensitivity to change in habitat. It ranges from a minimum of 0 
to a maximum of 9. A very common, widespread species that is highly tolerant to human distur- 
bance scores 0. In contrast, a range-restricted, threatened and sensitive endemic species scores 9 
(Samway, 2008). 


To obtain cursory information on some of the threats that each of our study sites were facing we 
used definitions and hierarchical classification as described by Salafsky et al., (2008). A Standard 
Lexicon for Biodiversity Conservation: Unified Classifications of Threats and Actions was used to 
describe and to categorize threats and actions. This method was chosen to allow for a consistent 
framework for threats to ensure consistent identification of shared threats across protected areas 
(Table S5 in Supplementary file). All impact scores (negligible, low, medium and high) were based 
on estimates whereby we interviewed managers, ecologists and conservationists working within 
Matobo National Park as well as communities in surrounding areas (sensu Ervin, 2003). 
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Data analysis 


The Simpson Diversity Index (Morris et al., 2014; Simpson, 1949) was used to determine the diver- 
sity for the adults and nymphs per each sample. After testing for normality using the Kolmogorov- 
Smirnov normality test and the Shapiro test, all data were found to deviate significantly from a nor- 
mal distribution. The differences their seasonal variations across the one-year period were compared 
using nonparametric tests. 


The relationship between diversity indices and the physico-chemical parameters were determined 
using correlation analyses. To analyze the association of different families to the environmental vari- 
ables at the dams studied, we used a canonical component analysis (CCA) using the vegan package 
in R (Oksanen et al., 2010). 


The relative abundance was calculated for each microhabitat by dividing the number of individuals 
observed in each microhabitat by the total numbers of individuals observed in that family. Indicator 
value (IndVal) analysis (Dufréne & Legendre, 1997) was used to identify if any of the species were 
especially associated with any of the microhabitats. The IndVal indicator attains the maximum value 
(1.00) when all the individuals ofa species are found in a single treatment (high specificity) and when 
the species occurs in all replicates of a treatment (high fidelity). The IndVal computations (‘labdsv’ 
package) were performed in R, version 3.5.0 (R Core Team, 2019). A randomization procedure as- 
sessed the significance of the indicator value of each species. Species that achieved a significant Ind- 
Val score were selected for interpretation. 


Results 
Nymph composition (Objective 1) 


Six Odonata families were identified from a total of 427 individual nymphs recorded during the en- 
tire sampling period. The highest abundance of Odonata (166 individuals) was recorded during the 
cold, dry season followed by the hot-dry season (138 individuals) and the least during the hot-wet 
season (123 individuals). Six families were recorded in the cold-dry season, whereas only five fami- 
lies were recorded in the hot-dry season and the hot-wet season. Based on the total number of indi- 
viduals recorded, Gomphidae was the most abundant family during the entire sampling period for 
the nymphs (166 individuals), followed by the Coenagrionidae and Libellulidae families with 128 and 
98 individuals recorded respectively. The Gomphidae family had the highest species abundance re- 
corded in the cold-dry season and the Hot-wet season, whereas the Coenagrionidae was the most 
abundant family during the hot-dry season (Table S1 in Supplementary file). 


Adult composition (Objective 1) 


A total of 35 species and seven families were identified from a total of 2401 individual adults re- 
corded, with the highest number of adults recorded during the hot-wet season (1294 individuals out 
of 35 species and seven families). However, six families and 33 species were found in the cold-dry 
season and, only six families and 28 species were recorded in the hot-dry season. The family Libel- 
ullidae had the highest abundance of adults recorded, followed by Coenagrionidae family during the 
entire sampling period. Brachythemis leucosticta recorded the highest number of individuals (315) 
during the entire sampling period followed by Trithemis arteriosa and Crocothemis erythraea with 
296 and 191 individuals, respectively. However, 7: arteriosa recorded the highest abundance during 
the cold dry season (85 individuals) and the hot dry season (84) whereas B. leucosticta recorded the 
highest abundance [177 individuals during the hot wet season]. 31 out of the 35 species identified 
recorded a low DBI score of 3 or less (Table S2 Table S3 & Table S4 in Supplementary file). 
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The results showed no significant difference in the diversity of nymphs [X2 (2, N = 45) = 2.9, p > 
0.05] among the three seasons during the entire study area, over the one-year period (Figure 2A). The 
results also showed no significant difference in the diversity of adults [X2 (2, N = 45) = 1.3, p> 0.05] 
among the three seasons during the entire study area, over the one-year period (Figure 2B). 
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Figure 2. Odonata species diversity among the three seasons sampled over a one-year period (May 2016 
to April 2017), for nymphs (A) and adults (B). Error bars are based on standard deviations. 
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Relationships between Odonata and abiotic factors (Objective 2) 


We found no significant correlation between diversity indices and species richness with the physico- 
chemical parameters (Table 1). Considering the relationship between odonate families and physico- 
chemical parameters (Figure 3), pH was positively correlated with Libellulidae and Platycnemididae, 
whereas total dissolved solids (TDS) and conductivity (EC) were positively correlated with Libellu- 
lidae, Platycnemididae and Macromiidae. On the other hand, the Gomphidae family was found to be 
negatively correlated to TDS, EC, and pH. We found a weak correlation between Coenagrionidae and 
water quality parameters. Conversely, we found no correlation between the family Aeshnidae and 


physicochemical parameters. 
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Figure 3. Canonical component analysis illustrating the associations between odonate nymphs and se- 
lect physico-chemical parameters collected from five dams. 
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Table 1: Spearman’s correlation coefficient between physicochemical factors and aquatic insect diver- 
sity and richness at Dams in the Matopos World Heritage site. 

















Parameter Taxa S_ | Simpson Index (1-D) 
Total Dissolved Salts —0.1303 | -0.1327 
pH 0.07736 | 0.14116 
Conductivity (us cm”) —0.1303 | -0.1327 














Gomphidae were found in higher percentages in sandy microhabitats (> 95%) and gravel (> 50%) 
than all other habitats sampled. Aeshnidae Families were found only where a detritus bottom sub- 
strate was present, with or without the presence of aquatic plants. Macromiidae and Platycnemididae 
were rare across all seasons (Figure 4). 
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Figure 4. Relative abundances of nymphs in five substrata during three sampling times in the Matobo 
National Park. 
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The indicator value analysis showed that Coenagrionidae, Libellulidae and Aeshnidae families were 
all significantly associated with detritus while the Gomphidae family was found to be significantly 
associated with sandy substrates (Table 2). 


Table 2: IndVal analysis of Odonata nymphs associated with specific habitat types. Only significant 
(p<0.05) habitat associations are presented. The significance of IndVal indices was assessed using 10,000 
Monte Carlo permutations. 

















Taxa Habitat | indval pvalue | freq 
Coenagrionidae Detritus | 0.22 0.001 53 
Libellulidae Detritus | 0.14 0.009 47 
Aeshnidae Detritus | 0.07 0.043 14 
Gomphidae Sand 0.49 0.001 56 























Threat assessments (Objective 3) 


Considering overall threats averaged across all dams, natural system modifications driven by fires 
and fire suppression as well as dams and water management, were identified as the most significant 
threats on Odonata around Lentic habitats in Matobo National Park [contributing ~ 18% of all as- 
sessed threats] (Table. S5 in Supplementary file). Other factors that showed a significant impact in- 
clude Agriculture expansion and intensification, specifically, the livestock from surrounding the 
communities which are left to drink water and graze around the dams within the National Park and, 
residential development (for Tourism and recreation areas), which contributed ~ 16% of all assessed 
threats. Habitat shifting and alterations, human intrusions and disturbances through recreational 
activities, pollution from garbage and solid waste and, invasive alien plant species (Figure S1 in 
Supplementary File) were also identified as the threats to Odonata around Lentic habitats in Matobo 
National Park, contributing ~ 10% of all assessed threats. It was also observed that over-exploitation 
of aquatic species is also a potential threat to Odonata around Lentic habitats in Matobo National 
Park contributing ~ 7.5% of all assessed threats (Figure 5, Table S5 in Supplementary file). 
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Figure 5. Direct threats impact of odonates at lacustrine systems in Matobo National Park. Only impacts 
> 5% are present. 
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Discussion 


We investigated the community assembly of Odonata species occurring in lacustrine systems in a 
World Heritage Site (Matobo National Park). Our results revealed that the nymphs and adults of the 
family Libellulidae and Coenagrionidae were the most abundant, and over 80% of the species re- 
corded had low DBI score indicating that the Odonata assemblage in lacustrine systems are mostly 
composed of wide spread generalist species, as documented in parts of the African continent (Ma- 
fuwe & Moyo, 2020). Additionally, Agricultural expansion and intensification, Residential and Com- 
mercial development, Habitat shifting and alteration through climate change and, human encroach- 
ment and natural system modifications were identified as major threats to odonate and other 
freshwater organisms in Matobo National Park. Contrary to our predictions, we found no significant 
differences in the diversity of nymphs and adults across all seasons. Our study makes a huge contri- 
bution to Odonata research and conservation initiatives in that it provides baseline data that can be 
used to monitor Odonata populations and distributions in aquatic habitats in the Matobo world heri- 
tage site. 


Nymph and Adult compositions 


Gomphidae were the most abundant nymphs, while Libellulidae were the most abundant adults. 
Specifically, the species, Brachythemis leucosticta in the Libellulidae family recorded the highest 
number of individuals during the entire sampling period followed by Trithemis arteriosa and Croco- 
themis erythraea. The abundance of the species mentioned above is tenable, considering that these 
three species are widespread in Africa (IUCN, 2020; Dijkstra & Clausnitzer 2014; Damm & Hadrys 
2012). The results also revealed that over 80% of the odonate assemblages at the dams in Matobo 
World heritage site have a low DBI score of 3 or less. This also indicates that the assemblage com- 
posed of mostly of common, widespread generalists odonate species which are tolerant to human 
disturbance (Simaika & Samways, 2009; Samways, 2008). On the other hand, some sensitive species, 
which are not very tolerant of disturbed areas and were last recorded in Matopos in the year 1980. 
e.g., Agriocnemis pinheyi (Pinhey wisp) and Bradinopyga cornuta (Horned rock dweller) were not 
recorded during this study. However, we only considered a short temporal scale of one year between 
2016 and 2017. 


Considering seasonal variation in diversity, there were no significant differences in the diversity of 
nymphs and adults across all seasons. The lack of seasonal variation may be explained by the fact 
that dams act as reservoirs for rainwater (Kosterin, 2010; Loureiro & Pontes, 2013) and they provide 
microhabitats, perch sites and artificial breeding sites that can support many species all year round. 
As such lentic habitats may support the same diversity of species across all seasons. However, with 
limited sampling the data in this study need to be interpreted with caution as sampling was done over 
one year between May 2016 and April 2017. It was reported that Zimbabwe experienced a severe 
drought in the year 2016, which was declared a national disaster with the drought being so severe in 
the southern parts of Zimbabwe (Government of Zimbabwe, 2017). This was also followed by some 
heavy flooding from the effects of cyclone Dineo in early 2017. During alternation of these extreme 
climates, odonates are known to sometimes alter their behaviour. For instance, some species may be 
present during wet years, but completely vanish where there have been ecosystem disruptions and 
extreme disturbances, only to appear and recover again when seasons have stabilized. Most of the 
freshwater habitats are also degraded or may also disappear during alternation of these extreme 
events (Abrahams et al. 2013) and recovery of odonate populations where habitats have been de- 
stroyed is very slow (Samways, 2010). Frequent sampling, over a longer period may therefore be able 
to reveal other patterns that were not revealed during this study. 
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Relationships between Odonata and abiotic factors 


Diversity and species richness were not correlated to physicochemical parameters. These findings 
are similar with results in a study by Ab Hamid, Salmah, & Nurul Huda (2016), who found no sig- 
nificant correlation between physico-chemical parameters (e.g. pH and TDS) and total abundance of 
Odonata nymphs. Notwithstanding, we only collected three abiotic variables thus, we possibly 
missed some of the variables that may have been critical as drivers of odonate nymph diversity. For 
example, dissolved oxygen, air temperatures, water temperature and flow rate have been documented 
as factors that significantly impact odonate diversity (Adu & Oyeniyi, 2019). Similarly, our taxo- 
nomic resolution for the odonate nymphs may have impacted the relationships between physico- 
chemical parameters and odonate nymphs. Specifically, we identified the odonates nymphs to mor- 
phospecies. Nevertheless, Odonata nymphs are difficult to determine up to species level, especially 
at the early instar level. This is because the characteristics used to identify nymphs are not fully de- 
veloped for early instars (Tennessen, 2019). Fraser (1956) states that the best method is to rear them 
in for correct identification up to species. 


We found the physicochemical quality parameters were related to odonate families. For example, pH 
was positively correlated with Libellulidae and Platycnemididae, whilst total dissolved solids (TDS) 
and conductivity (EC) were positively correlated with Libellulidae, Platycnemididae and Macromi- 
idae. Many species of Libellulidae and Platycnemididae are known to be salt-tolerant, and many are 
generalists that can inhabit in almost all kinds of water, whether acidic, alkaline, brackish or saline 
(Lakhiar & Panhwar, 2015; Saeed et al., 2016). However, the Gomphidae family was found to be 
negatively correlated to TDS, EC and pH. According to Ab Hamid, Salmah, & Nurul Huda (2016), 
gomphids thrive better in freshwaters compared to brackish waters. Additionally, Gomphidae are 
sensitive to pollution and environmental disturbance (Bouchard, 2004). 


The results showed the highest percentage of Gomphidae nymphs were found in the sandy substrates 
followed by the gravel bottom substrates across all seasons (Figure 4). The Indval analysis also 
showed the Gomphidae family to be significantly associated with sandy bottom substrates. Gomphi- 
dae as sedentary burrowers are found where there is sand, silt, mud or gravel to burrow into (Leipelt 
et al., 2010; Palatov, Chertoprud, & Frolov, 2016; Thorp & Covich, 2015). 


The Coenagrionidae and Platycnemididae nymphs were recorded in aquatic plants. Coenagrionidae 
and Platycnemididae species are known to be found amongst submerged vegetation, such as willow 
roots and reeds (Thorp & Covich, 2015). Females lay their eggs in plant tissue, and the nymphs cling 
on to the submerged plants (Thorp & Covich, 2015). Coenagrionidae nymphs were also found in 
detritus substrates possibly because Coenagrionidae nymphs are known to survive in detritus even 
when water sources dry up and some dams like Maleme dam were observed to be very sluggish with 
a lot of detritus especially during the dry season when they were drying up. 


Aeshnidae was only found where there was a detritus bottom (with or without aquatic plants present). 
Aeshnids are known to prefer a wide range of habitats, including brackish and sluggish waters char- 
acterized by mud and detritus (Hawking et al., 2004). Aeshindae are also known as climbers adapted 
to living on vascular hydrophytes or detrital debris (e.g., overhanging branches, roots and vegetation 
along streams, submerged brush in lakes) with modifications for moving vertically on stem-type 
surfaces. In this study, a large number of samples were collected from detritus bottom substrate. The 
muddy substrates may be a result of the deposits of animal waste from the livestock and wildlife 
drinking from the dams. The organic detritus was also composed of fine silt mud possibly brought 
down by rivers and land runoff. 


Macromiidae was only recorded during the cold, dry season where there was mud and detritus. Mac- 
romiidae are known to occur in mud, silt, gravel, rocks and litter in the benthic zone (Laughlin et al., 
2018), however, only a few species of Macromiidae species are found in Africa and most reproduce 
in flowing waters, which can explain why only one individual was recorded at these lentic habitats 
throughout the study period. 


Libellulidae was found across most of the microhabitats across all seasons. Many species in the Li- 
bellulidae are widespread, and the family has many generalist species which are known to occupy a 
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wide range of still or sluggish waters; including streams that dry to pools in the outback, riverine 
pools, lagoons, lakes, boggy seepages, swamps and permanent ponds (Samway 2008; Dijkstra & 
Clausnitzer, 2014). 


Threat assessments and Conservation of Odonates 


Considering threats to Odonata around lentic habitats in Matopos National Park we found natural 
system modifications to be the most significant threats. The likely modifications in MNP are likely 
related to hydroperiods (water level fluctuations) of dams that are affected by rainfall patterns in 
Southern Africa, as documented by other researchers (sensu Deacon, Samways, & Pryke, 2018; 
Briggs et al. 2019; Jooste, Samways, & Deacon, 2020). According to IUCN (2020), the ‘natural sys- 
tem modifications’ is a threat category which covers a wide range of activities that convert or degrade 
habitats largely as a result of human management. The dominant threat in this category for aquatic 
systems is damming, which fragment river systems, alter flows and sedimentation patterns, and 
cause mortality and decline in aquatic organisms (IUCN, 2020). For example, large hydroelectric and 
flood control dams which flooded 102,000 hectares of wetlands in the Columbia Basin possibly 
eliminated populations of almost all dragonfly species on the Basin list and; the modification of the 
shores of the Kootenay and Columbia rivers through dam construction has eliminated many habitats 
for Argia emma, Ophiogomphus occidentis and Macromia magnifica. Similar natural system modi- 
fications have been implicated in the decline of other freshwater invertebrates at many locales (Vin- 
son, 2001). 


Agriculture expansion & intensification, specifically, the livestock from surrounding the commu- 
nities which are often left to drink water and graze around the dams within the National Park (K. 
Mafuwe, Personal observation, November 2016), and development of tourism facilities and recre- 
ation areas, also contributed approximately 16% of all assessed threats. Livestock grazing in the 
national park has been observed and mentioned by a number of tourists and visitors to the park as 
early as the year 2012 (Newsday - September 4, 2012; Zimfieldguide.com, n. d). Cattle grazing and 
trampling around aquatic habitats, has been known to increase habitat alteration by destroying 
vegetation, reducing plant cover, destroying soil structure and, increasing surface runoff of water 
and siltation (Dunne, Western, & Dietrich, 2011; Pelinson, Garey, & Rossa-Feres, 2016) e.g. ob- 
served around Mesilume dam (Figure S2 in Supplementary File). Dams around Matopos National 
Park are used as camping locations and for recreational activities such as picnic and fishing sites by 
tourists. Direct and indirect human activities including tourism and recreational activities are 
known to cause destruction of the freshwater habitats through pollution, littering and destruction 
of vegetation. This consequently has a significant impact on aquatic organisms. Similar impacts 
may also result from the human intrusions and disturbance and, pollution threats which were ob- 
served to contribute about 10% of all assessed threats. For example, Spyra & Strzelec (2019) con- 
firmed that recreational activities alter some characteristic of environments such as altering the 
substrates for perching such as leaf deposits and plants, impacting on benthic invertebrate popula- 
tion. 


Habitat shifting & alteration and invasive alien species were also identified as the threats to Odonata 
around Lentic habitats in Matobo National Park and also contributed ~ 10% of all assessed threats. 
Similarly, invasive alien species have been recognized as a threat to Red-Listed dragonflies in coun- 
tries like South Africa’s (Samways & Taylor, 2004). Elsewhere, Garcia et al. (2010) also notes that 
invasive species threaten freshwater species. Further, anthropogenically-induced climate change 
increases the frequency and occurrence of droughts which subsequently result in habitat shifting & 
alterations. The fluctuating water levels not only affect the quality of the water for example by chang- 
ing the acidity levels but, also impact on the quality of the surrounding vegetation, possibly resulting 
in a loss of some Odonata species, e.g. Lestes dryas Kirby, 1890 a species of dragonfly which de- 
clined due to climate change which had an impact on many habitat factors, in England and Ireland 
around the year 1980 (Moore, 1980). The fluctuations in temperature may also influence changes in 
geographical distributions, through a shifting of species’ fundamental niches, changes in the distri- 
bution of suitable habitat and variation in the dispersal ability of species (Hassall & Thompson, 
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2008). We also recorded that over-exploitation of species as a potential threat to Odonata around 
Lentic habitats in Matobo National Park contributing ~ 7.5% of all assessed threats. More specifi- 
cally, the use of gillnets for fishing nets was observed at Maleme dams during the study period. The 
use of fishing nets is particularly destructive to the none target species and their freshwater habitats 
and therefore the impact extends further than just to the fish. Considering the threats to the Matobo 
National Park, we recommend that the Zimbabwe National Parks and Wildlife Management Author- 
ity embarks on more stringent regulations on fishing practices, invasive species management, and 
abstraction of natural resources in the Matobo Nation Park. Future construction of tourism facilities 
should also aim at minimizing the potential impacts of land-clearing and modifications of the land- 
scape during the construction of tourism facilities, such as removal of vegetation, changes in water 
flows, and destruction or alteration of breeding and feeding grounds of freshwater species. In addi- 
tion, such habitat alterations could to be compensated by additional habitats. Even further, wetlands 
could be promoted through tourism facilities. 


Summarily, this study has documented the seasonal diversity and distribution of odonates (nymphs 
and adults) in an understudied protected area. A number of species and families recorded in this 
study are known to have wide range and very tolerant of polluted and disturbed waters e.g., Tri- 
themis arteriosa and Crocothemis erythraea and many other species of the Libellulidae and Plat- 
ycnemididae families. Additionally, over 80% of the species recorded have a low DBI score indi- 
cating that the Odonata assemblage at the lacustrine systems is mostly comped of widespread 
generalist species. Some sensitive species, which are not very tolerant of disturbed areas and were 
last recorded in Matopos in the year 1980 were also not recorded during this study. However, only 
a short temporal scale (2016/2017) was considered. Considering Zimbabwe experienced severe 
weather extremes during the study period, long monitoring data over long temporal scales are 
warranted. The threats on Odonata species and the freshwater habitats, identified in this study 
provide a basis for call to actions to protect freshwater biodiversity in Zimbabwe. Because Odo- 
nata species have been identified as key indicator species, habitat protection and restoration are 
the most important steps towards conservation of these species. It is necessary to maintain the 
structural integrity of both larval and adult habitats, i.e. water bodies and their surrounding land- 
scape. 


Broadly, considering the environmental and anthropogenic pressures in water resources and the 
world over, our findings stress the need for more focused baseline studies documenting diversity and 
abundances of freshwater species; these baseline studies will allow ecologist and environmental 
managers to predict the effect of future environmental and anthropogenic perturbations. Further, 
considering that odonates are not listed in many conservation strategies in most African countries, 
there is need for baselines studies to explore the conservation strategies for small taxa animals such 
as odonates. Our results will be useful to conservation managers in the Matobo National Park and 
other protected areas in Sub- Saharan Africa and the world. 
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